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In this work, we show that the zero field excitation spectra in the quantum spin ice candidate 
pyrochlore compound Yl^T^Or is a continuum characterized by a very broad and almost flat 
dynamical response which extends up to 1 — 1.5 meV, coexisting or not with a quasi-elastic response 
depending on the wave-vector. The spectra do not evolve between 50 mK and 2 K, indicating that 
the spin dynamics is only little affected by the temperature in both the short-range correlated and 
ordered regimes. Although classical spin dynamics simulations qualitatively capture some of the 
experimental observations, we show that they fail to reproduce this broad continuum. In particular, 
the simulations predict an energy scale twice smaller than the experimental observations. This 
analysis is based on a careful determination of the exchange couplings, able to reproduce both the 
zero field diffuse scattering and the spin wave spectrum rising in the field polarized state. According 
to this analysis, Yb 2 Ti 2 07 lies at the border between a ferro and an antiferromagnetic phase. 

These results suggest that the unconventional ground state of Yb 2 Ti 207 is governed by strong 
quantum fluctuations arising from the competition between those phases. The observed spectra 
may correspond to a continuum of deconfined spinons as expected in quantum spin liquids. 


I. INTRODUCTION 

Understanding, characterizing and classifying novel 
states of matter is one of the main goal of the research in 
solid state physics. In particular, systems where the ther¬ 
mal or quantum fluctuations are able to melt long range 
order, the so-called spin liquids, draw a lot of attention 
since they generally go beyond Neel paradignP. 

The pyrochlore lattice, made of corner-sharing tetra- 
hedra is the archetype of a three-dimensional frustrated 
lattice and has proven during the last years to be a rich 
playground for studying such spin liquid stated. Among 
the variety of possible pyrochlore systems, the compound 
Yb2Ti207 has been presented as one of the possible real¬ 
izations of a quantum spin liquid ground-state. More pre¬ 
cisely, it has been proposed as a candidate for the quan¬ 
tum variant of the spin ice state observed in Hu^TUOr 
and Dy2Ti207. In these classical spin ices, the Ising-like 
anisotropy of the magnetic moments along the local (111) 
directions, combined with an effective ferromagnetic in¬ 
teraction induce a macroscopically degenerated ground 
state characterized by the local ice-rule^ (in each tetra¬ 
hedron, two spins point in and two spins point out). In 
Yb2Ti207, the Yb 3+ magnetic moment shows a weak 
XXZ planar anisotropy perpendicular to the local (111) 
direction^®. However, the strongly anisotropic interac¬ 
tions tensor^, whose main component is ferromagnetic, 
induces a local constraint analog to the ice-rule. 

The static and dynamical magnetic correlations in 
Yb2Ti207 have been investigated by means of neutron 
scattering experiments. Short range correlations settle 
around To ~ 2 K, giving rise to rods of diffuse scattering 
along high symmetry directions 9 11 . At lower tempera¬ 


ture, a phase transition is observed in specific heat mea¬ 
surements around 0.2 m but the critical temperature 
depends on the nature of the sam ple s (polycrystal or sin¬ 
gle crystal) and their quality 11 13 14 . This transition was 
shown to be first order 14 "Uj in both single crystals and 
powder samples but its nature remains debated. Several 
studies (including neutron scattering and muon spin re¬ 
laxation measurements) evidence the stabilization of an 
ordered ferromagnetic state, with a reduce d stat ic mag¬ 
netic momentP^KEl, while others do not ! 10 * 14 * 1 -^^. 

To understand the underlying microscopic mechanism 
leading to such a non-conventional ground-state, several 
studies attempted to determine the exchange couplings, 
trying to reproduce the experimental results. Different 
sets of interactions were obtained, depending on the fit¬ 
ted quantity: diffuse scattering 9 16 or field induced spin 
wave excitation^. The latter parameters allow repro¬ 
ducing specific heatPS and were further refined using 
terahertz spectroscopy measurements 22 . These parame¬ 
ters place Yb2Ti2 07 in a ”splayed ferromagnetic” phase, 
quite far from the ” quantum spin liquid” phase of the 
theoretical phase diagram 23 2 -4 

Since the precise nature of the elementary excitations 
generally depends on the ground-state itself, we propose 
in this article to address the unconventional magnetic 
properties of Yb2Ti207 through a detailed study of the 
spin dynamics on a single crystal in the zero field cor¬ 
related phase below 2 K, which has been little explored 
up to now 11 . In section |n| we first present macroscopic 
measurements performed on a piece of our single crys¬ 
tal, which show evidence for a transition at 175 mK. By 
means of inelastic neutron scattering, we show that the 
excitations spectrum is a continuum characterized by a 
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very broad and almost flat dynamical response extending 
up to 1 — 1.5 meV, and coexisting or not with a quasi¬ 
elastic response depending on the wave-vector. These 
excitations do not depend on the temperature below the 
stabilization of short-range correlations (i.e. 2 K), and, in 
particular, entering the ordered phase does not affect the 
spectra. Then, in section m we compare these inelastic 
data with classical calculations combining Monte Carlo 
and spin dynamics simulations. We show that, although 
the qualitative features of the spectra are reproduced, 
calculations predict an energy scale twice smaller than 
the experimental observations. These calculations are 
based on our determination of a reliable and robust set 
of exchange couplings, able to reproduce both the spin- 
wave spectrum of the field polarized phase as well as the 
diffuse scattering in zero field. The obtained couplings do 
place Yb 2 Ti 2 07 in a ferromagnetic phase, yet very close 
to the boundary with an antiferromagnetic phase. Then, 
although Yb 2 Ti 2 07 is not a canonical quantum spin liq¬ 
uid as initially proposed, we propose that the proximity 
of competing ferromagnetic and antiferromagnetic corre¬ 
lations restore strong quantum fluctuations— and is at 
the origin of the anomalous static and dynamical behav¬ 
iors observed in our experiments. 


II. EXPERIMENTAL 

Our measurements were performed on a single crystal. 
First, a polycristalline sample of Yb 2 Ti 2 07 stoichiometry 
was synthesized from Yb 2 03 and Ti0 2 starting powders 
by solid state reaction. To obtain the pyrochlore phase 
several thermal treatments at 1400 °C with intermediate 
grindings were necessary. The progress of the synthe¬ 
sis reaction was followed by powder X-ray diffraction. 
For the last thermal treatment the powder was shaped 
as a cylinder of 5 mm diameter and 90 mm length and 
the obtained road was used as feed road for the single 
crystal synthesis. Crystal growth was performed using 
the floating zone technique in a four-mirror optical im¬ 
age furnace NEC SC-N15HD. The obtained crystal was 
then annealed under 0 2 gas flow for two days. 


A. Thermodynamic measurements: Magnetization 
and specific heat 

To characterize the macroscopic properties of our 
Yb 2 Ti 2 07 sample, we have performed thermodynamic 
measurements, magnetization and specific heat, on a flat 
disk sample (mass 244 mg), cut from the single crystal 
used for neutron scattering experiments. These measure¬ 
ments were performed at the Institut Neel in purpose 
built experiments equipped with 3 He- 4 He dilution refrig¬ 
erators. Magnetization was measured in a SQUID mag¬ 
netometer as a function of field and temperature down to 
90 mK, together with ac susceptibility measurements 2 -^. 
Specific heat was measured by the relaxation method 



FIG. 1. (Color online) Specific heat C (red squares) and mag¬ 
netization M/H (blue dots) vs temperature T. The specific 
heat is the long time specific heat (see text) and the magne¬ 
tization was measured in a field of 3.83 mT, with a step of 10 
mK every 1350 s. 


down to 70 mK in several applied fields up to 100 mT. 
The high sensitivity and fast response of the thermome¬ 
ters used in this experimental set-up as well as their high 
stability allow to measure both short and long time heat 
relaxation 27 . The magnetic field was applied in an arbi¬ 
trary direction, parallel to the plane of the disk. 

The magnetization exhibits a Curie-Weiss temperature 
of 0.5 K, and reaches 1.65 /ig/Yb at 90 mK an d 3 T w hich 
is consistent with previous measurement j511 S l 18 l 28 l 29 J The 
specific heat presents a broad maximum around 2.5 K 
characteristic of short-range order correlations in agree¬ 
ment with previous results 12 . When decreasing the 
temperature, these measurements show a transition at 
Tc = 175 mK, in both specific heat and magnetization, 
which is evidenced by a peak in the specific heat, an 
upturn in the magnetization (see Figure [l]), and the on¬ 
set of the out-of-phase ac susceptibility (not shown). As 
previously observed^, the magnetization shows a small 
hysteresis between the cooling and warming curves be¬ 
low the transition, indicative of a first order transition. 
Nevertheless, although the peak in specific heat is rather 
narrow, the increase of the magnetization when decreas¬ 
ing the temperature below the transition is smooth, and 
fails to saturate, even at 90 mK. This result suggests that 
the magnetic transition expands on a broad temperature 
range and that no net spontaneous ferromagnetic mo¬ 
ment is stabilized in this sample down to 90 mK. 

When a small field is applied {fi^H < 100 mT), the 
transition persists and is shifted to higher temperatures, 
confirming the ferromagnetic nature of the transition (see 
Figure [2] for specific heat). It is worth noting that below 
1 K, the heat relaxation is not exponential, resulting in 
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FIG. 2. (Color online) Specific heat C vs temperature T for 
several magnetic fields. Open symbols correspond to the long¬ 
time specific heat (t ~ 600 s) and full symbols to the short 
time specific heat (fw 9 s). 


two contributions for the specific heat: the short time one 
(characteristic time of about 9 s) and the long time one 
(characteristic time of about 600 s). These contributions 
show qualitatively the same temperature dependence in 
zero magnetic field, which evidences that dynamics with 
different time scales exist in the system above and be¬ 
low the transition. When increasing the magnetic field, 
the peak in the short time response is suppressed more 
drastically than the long one (See Figure [2]). 


B. Neutron scattering 

The neutron scattering measurements were performed 
on a large Yb 2 Ti 207 single crystal grown with the float¬ 
ing zone technique, as detailed above. The crystal was 
aligned in the (hhO) — (00^) scattering plane and cooled 
down to 50 mK in a dilution fridge. Note that due to the 
large size of the crystal, its poor thermal conductivity 
at very low temperature, and the heat provided by the 
neutron beam, the true temperature of the crystal might 
be slightly larger. Special care to the thermalisation of 
the sample was taken during the experiments: a wait¬ 
ing time of a few hours was used between the measure¬ 
ments to ensure a stable temperature inside the sample. 
Polarized and unpolarized experiments were respectively 
conducted on the 4F1 and 4F2 triple-axis spectrometers 
(LLB, France), with final neutron wave-vectors kf =1.3 
A -1 , resp. kf =1.15 A -1 , yielding an energy resolu¬ 
tion of A r = 150 /ieV, resp. A r = 85 /ieV. A cooled 
beryllium filter was placed in the incident beam to re¬ 
move high order contaminations. The polarization anal¬ 
ysis allows measuring the spin-spin correlation functions 
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FIG. 3. (Color online) (a) Elastic intensity map in the ( h h i) 
scattering plane at T = 0.05 K. White lines correspond to 
the Brillouin zones boundaries. Black dashed lines denote 
the directions probed in the inelastic neutron scattering ex¬ 
periments (see section [il B 2| ). (b) Elastic intensity (full dots) 
along the (h h 2 — h) direction obtained by integrating Figure 
(a) over 5Q = 0.4 rlu in the perpendicular direction. The lines 
with open symbols show calculated elastic scattering functions 
obtained from classical spin dynamics simulations for differ¬ 
ent parameter sets J 2 = —0.26, —0.29, —0.3 1 meV (gray), 
—0.326 meV (blue) and T = 0.4 K (see section HI E| . (c) Spin 
flip elastic intensities (full symbols) along the (h h 2 — h) di¬ 
rection at T = 0.9 K with polarization applied along y (blue 
dots) and z (red squares). The blue and red lines with open 
symbols correspond to calculated elastic scattering functions 
Sy (Q,u; = 0) and S z (Q,w — 0) obtained from classic al spin 
dynamics simulations for J 2 = —0.32 meV (see section HIE). 


S„(Q , lS) and S z (Q,cj) where the subscripts y and z in¬ 
dicate that the spin components are perpendicular to Q 
within the scattering plane and along the vertical axis 
respectively. Both are measured in the spin flip chan¬ 
nel, with polarization (direction of the incident neutron’s 
spin) P applied along z and y respectively. The flipping 
ratio was about FR=42 at high temperature but slightly 
decreased below T = IK, down to FR=33 at T = 0.65 K. 
Because of technical difficulties during the polarized neu¬ 
tron experiment, the lowest reachable temperature was 
T ~ 650 mK. Note that the magnetic intensity measured 
in unpolarized experiments is the usual spin-spin corre¬ 
lation function S(Q,u) = S^Q,^) + S z (Q,w). 


1. Static properties 

Elastic data collected at T = 0.05 K are shown in Fig¬ 
ure [3j revealing the same qualitative features as the en¬ 
ergy integrated diffuse scattering described at len gth in 
prior works above the critical temperaturd^ 11 16 20 . It is 
characterized by a line of scattering along {h h h) accom¬ 
panied by another branch from Q = (1,1,1) to (2,2,0) 
ending with a large spot at Q = (2, 2, 0). The presence of 
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spectral weight around both (1,1,1) and (2, 2,0) is better 
evidenced in Figure [3] (b), where the elastic intensity is 
plotted versus wave-vector along the (h h 2 — h) direction. 
While these rods of diffuse scattering are undoubtedly 
related to the presence of strongly anisotropic exchange 
interaction^- 9 , their origin is still unclear. First ascribed 
to a dimensionality reduction of the static correlations 11 , 
the ground-state has more recently been shown to be fully 
3D in the region of parameter space where Yb 2 Ti 2 07 is 
supposed to ikP. 

The spectral weight along the rods is particularly in¬ 
tense around the (2, 2, 2) Brillouin zone center, as ex¬ 
pected in the presence of ferromagnetic correlations : in 
the limit of a long range collinear ferromagnetic order, 
the strongest intensity is obtained for (2n, 2n, 2 n) Bragg 
peaks, while the (2 n + 1, 2n + 1, 2n + 1) ones are around 
three times weaker and the (2n, 2n, 0) peaks are extinct. 
A finite spectral weight at Q = (2, 2, 0) is characteris¬ 
tic of antiferromagnetic spin arrangements, which points 
out the coexistence of both ferro- and antiferromagnetic 
correlations in the system. 

These results reveal that the short range correlations 
still remain at very low temperature and coexist with the 
long range ferromagnetic order. This might be related to 
the fact that, in our sample, the magnetization increases 
smoothly below the transition temperature. 

Finally Figure[3](c) displays data in the same direction 
measured using polarized neutrons at higher temperature 
(T = 0.9 K). It points out that the signal around Q = 
(2, 2, 0) is polarized along z, while the one around (1,1,1) 
has similar intensity in both the y and z channels. The 
present data are thus consistent with the one obtained in 
Ref. 46, 


2. Dynamic properties 

Series of inelastic data were then collected along sev¬ 
eral high symmetry directions (h h 2), (0 0 f), (1 1 Q, 
(h h 2 — h) and (h h 3 — h) at the base temperature of 
T = 0.05 K. For all measured wave-vectors Q (see Figure 
[3ja), black dots), the spectra show above the elastic line 
a very broad signal extending up to about 1 — 1.5 meV 
(see Figure [4]). Despite the rather good energy resolu¬ 
tion A r =85 /ieV, no well-defined collective excitations 
could be observed. This is in sharp contrast with the 
resolution-limited spin- wave s observed in the field polar¬ 
ized phase of Yb 2 Ti 2 O 7 ^ 10 l 

Two qualitatively different behaviors can be observed 
depending on the wave-vector position relative to the 
rods of scattering. The spectra taken at Q vectors lo¬ 
cated on these rods are dominated by quasi-elastic fluc¬ 
tuations, whereas away from them, the signal appears 
inelastic. This feature is evidenced on the raw data 
presented in Figure |4j displaying directions (h h 2), 
(h h 2 — h) and (0 0 i) on the left, middle and right 
panels respectively. On the left panel, while ” flat-toped” 
at Q = (0,0,2), the signal becomes quasi-elastic on ap- 
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FIG. 4. (Color online) Excitation spectra along high symme¬ 
try directions (h h 2), (0 0 £), and (h h 2 — h) at T = 0.05 K. 
The blue area stands for the background. Dashed red and 
blue plain lines are the result of the fit considering a dual re¬ 
sponse consisting in quasi-elastic and inelastic contributions. 
Each contribution is represented in red while the total is rep¬ 
resented in blue. 


proaching the rod position at (2, 2, 2). The same feature 
is observed for the scans taken along (h h 2 — h) : the 
line-shape of the spectra becomes quasi-elastic-like while 
going through the rod position at Q = (1,1,1) and at 
(2,2,0). Finally, along (0,0 ,£), i.e. away from the rods, 
an inelastic line-shape is observed. 

A phenomenological fit of the data can be performed 
considering a dual response consisting in quasi-elastic 
and inelastic contributions multiplied by the detailed bal¬ 
ance factor (see Figure [4j plain and dashed lines). Since 
the experimental resolution is one order of magnitude 
smaller than the observed excitations, the fitting function 
was not convolved with the experimental resolution func¬ 
tion. The fitting procedure is as follows. First, the width 
Yqe — 0.3 meV of the quasi-elastic response is evaluated 
using a Lorentzian profile at different wave-vectors along 
or close to the rods of scattering. Yqe is then considered 
as wave-vector independent. In addition, a broad Gaus¬ 
sian peak is used to model the inelastic response. The fit 
is able to converge only if the width of this inelastic peak 
is fixed to a given value. The free parameters of the fit 
are thus the intensity of both quasi-elastic and inelastic 
contributions as well as the peak position. This modeling 
shows qualitatively the different behaviors close to and 
away from the rods of scattering, respectively dominated 
by quasi-elastic and inelastic response. Note that, in such 
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an analysis, the inelastic contribution is peaked at about 
0.5 meV and seems not to disperse. 

These inelastic data, and especially the flat energy de¬ 
pendence, are for the most part consistent with earlier 
inelastic results obtained below Two = 0.7 meV around 
two positions only, Q ~ (1.75,1.75, 0.5) and (1.5,1.5,1.5) 
(see Figure 4 in Reference mj). The present results 
demonstrate that the on- and off-rods behaviors can be 
generalized throughout the Brillouin zone. Furthermore 
they show that the dynamical response extends up to 
1 — 1.5 meV (see Appendix |A| for a further detailed anal¬ 
ysis). 


3. Temperature dependence 

Interestingly, the spectra recorded at higher tempera¬ 
ture do not show any change compared to the low temper¬ 
ature data, indicating that the spin dynamics is only little 
affected by the increase of temperature up to 850 mK. 
This remains valid whether Q lies close to the rods of 
diffuse scattering or not (see Figure [5] (a,b,d) and (c,e) 
respectively). Although the lack of temperature depen¬ 
dence of the “off-rod” scattering is consistent with Ref. 
QU the reported “on-rod” depletion at energies less than 
0.2 meV is not observed in our data. The spectra remain 
quasi-elastic. As a result we cannot conclude on a slow¬ 
ing down of the spin fluctuations on entering the ordered 
phase. This difference may be explained by the low tran¬ 
sition temperature observed in our sample, as denoted in 
section III A1 

Inelastic polarized neutron measurements have also 
been performed at higher temperature (T = 2 and 
4.5 K) at Q = (1,1,2). The data, shown in Fig¬ 
ure [6j first confirm the magnetic nature of the inelas¬ 
tic signal. At T = 4.5 K, the full polarization anal¬ 
ysis points out that the S(Q,uj) y contribution is a lit¬ 
tle bit more intense than the S(Q,co) z . Surprisingly, it 
also appears that S(Q,uj) z (red curve) vanishes above 
h^max — 0-75 meV while S(Q,uj) y (blue curve) extends 
up to ~ 1.1 meV. A depletion of the spectrum is also ob¬ 
served at low energies, making the generalized suscepti¬ 
bility x"(Q,u;) = S(Q,cj)/ (1 + n(uj)) quite narrow, with 
a maximum around Tiuu = 0.6 meV (see Figure m ( c ))- 
Interestingly, this approximately well defined excitation 
observed at T = 4.5 K spreads out and becomes less de¬ 
fined with decreasing temperature (cf. Figure |6jb)). At 
T = 0.05 and 2 K, x"(Q,u;) is superimposed, showing 
that the spectrum does not evolve anymore below 2 K 
(see panel (c)). 


4- A possible continuum of excitations 

In summary, this study of Yb 2 Ti 2 07 points out very 
unconventional spin dynamics. Starting from the high 
temperature, the spin excitations spectrum broadens be¬ 
low Tq ~ 2 — 4 K. This is largely unexpected, since, in 
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FIG. 5. (Color online) Excitations spectra at several wavevec- 
tors close to (a,b,d) or away from (c,e) the rods of diffuse 
scattering, and temperature from T = 50 to 850 mK. 


general, the width of inelastic spectra, related to the rate 
of the spin fluctuations, tends to increase with increas¬ 
ing temperature. The crossover temperature To between 
the two regimes coincides with the specific heat bump 
observed around T = 2.5 kP, below which short range 
spin ice like correlations establish 21 and the rods of dif¬ 
fuse scattering rise 9 11 20 . At low temperature, no well 
defined excitations are observed but rather a continuum 
of excitations. Finally, entering the ordered phase does 
not affect the spectra. 


III. NUMERICAL SUPPORT 

These experimental results are now confronted to nu¬ 
merical calculations. First, a model of exchange interac¬ 
tions at play in Yb 2 Ti 2 07 is proposed. Based on RPA 
numerical simulation d 30 ! 31 !, we show that a large range of 
parameters accounts for the available neutron data ob¬ 
tained in the conventional polarized phase induced by 
applying magnetic fields of 2 and 5T. Reproducing the 
strong quasi-elastic scattering around Q = ( 2, 2,0) with 
Monte Carlo Spin Dynamics simulations 32 3 ^ ] allows how¬ 
ever to constrain further the set of parameters and deter¬ 
mine an optimal set of coupling constants. In zero field, 
these simulations show that Yb 2 Ti 2 07 is a canted ferro- 
magnet, which yet lies very close to an antiferromagnetic 
phase. Finally, RPA and classical spin dynamics calcula¬ 
tions are presented. Although the latter reproduce some 
of the experimental features, they both fail to explain the 
very large energy range of the experimental data. This 
points out the role of quantum fluctuations that may be 
amplified by the proximity to a competing antiferromag- 
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FIG. 6. (Color online) (a) Spin polarized neutron intensity 
S(Q,u) taken at Q = (1,1,2) and T = 4.5 K. Red squares, 
blue circles and green triangles respectively denote the spin 
flip contribution in the z, y and x polarization channels. 
The inset shows the calculated spectra, see section |mE| at 
T — 0.9 K. (b) Spin flip intensity with polarization along 
x at 4.5 (red squares) and 2 K (blue circles), (c) shows 
the imaginary part of the magnetic susceptibility x^Q?^) 
at Q = (1,1,2) for 4.5 (red squares), 2 (blue circles) and 
0.05 K (green triangles). 


netic phase, and which are not captured using a classical 
approach. 


A. Models and Hamiltonians 

The starting point is the widely accepted Hamiltonian 
for pyrochlore systems 2 : 

H — Hcef + H eX c + H Z . (1) 

Hz = 9j J2i H ' is the Zeeman term, with H 
the applied magnetic field and Ji the magnetic mo- 
ment at site i. H C ef = Ei J2k,q= 1 B k^k,i is ttle 
crystal electric field (CEF) Hamiltonian, with 0\ the 
Stevens operators and B / scalar parameters 37 . Sev¬ 
eral estimations of the B\ parameters may be found 
in the literatur e! 4 * 5 * 38 * 39 *. In the present paper, we 
use the CEF parameters (T> 9 , H 9 , H 3 , Bq, Bjj, B%) = 

(536,7752,-2942,830,671,739) K obtained in ref. M 
These parameters lead to a doublet CEF ground-state, 
separated from the first excited level by an energy gap 
A E ~ 700 K. The magneto-crystalline anisotropy of 
the ground-state doublet is XXZ-like, with an easy 
plane perpendicular to the local (111) directions. The 
g z /g±_ — 0.504 ratio denotes a weak in-plane anisotropy 
with Lande factors g z = 2.06 and g± = 4.09. H exc = 
is a bilinear coupling Hamiltonian, where 

the interaction tensor JCij couples next neighbors mag¬ 
netic moments J at sites i and j. Note that H exc gath¬ 


ers all possible physical (e.g. exchange and dipolar cou¬ 
pling restricted to nearest neighbors), whether they are 
isotropic or not. By symmetry arguments, the 9 cou¬ 
pling constants of the 3x3 tensor JCij are reduced to 
only fouffl Here, we assume an exchange tensor which 
is diagonal in the (a, b, c) frame linked with a R-R bond 6 : 

JVJ, = Y, J i (*<*«• + ICbbfM + /C c cf.4) ^ 

H,v=x,y,z 

—JC^V^ X Jj) 

Since the energy gap between the CEF ground-state 
and the first excited levels are order of magnitudes larger 
than the exchange interactions and the Zeeman term, it 
is possible to define effective spin 1/2 operators, denoted 
by by projecting the full moment onto the CEF 
ground-state doublet. As a result, an effective Hamilto¬ 
nian can be defined in terms of new anisotropic couplings 
between the spin components: 

ff eff = X S 'T S ,’ (2) 

ij 

A popular convention used in RcffS. l23l 1241 and [40] con¬ 
sists in using (J±£, J z ±, J zz ) defined as: 


H e & = X s i - J ± ( s + s 7 + s,-st) 




+J±± (7^S+S++ 7 *.S-ST) 

+Jz± [S^ (GjS+ + QjSj ) + i++ j] 

where 7 ijXij are c-number J 8 I 23 I 24 [ Note that ” sanserif’ 
notations refer to local bases. (J±±, J±, J z ±, J zz ) are re¬ 
lated to JC a ,b,c ,4 by the following relations 

v2 JC a - 2 JC c - 4 /C 4 


Szz = K 


J± = -A 


2 2 JC a - 3 JC b -K c + 4 /C 4 


J z ± — Ax X z 

Aj_ 


12 

JC a + JC c — JC 4 
3y/2 

2 2 JC a + 3 JC b - JCc + 4/C4 

12 


with A^ 2 

- L ’ z gj 

In the following, we will follow the alternative con¬ 
vention of Ref [ 8 ] and (25, using the set of effective pa¬ 
rameters Ji, 2 , 3,4 (correspondence between the different 


conventions are 

given in Appendix |bJ) : 

Ji = 

h 

2J±± + 4J± + 2\/2J z ± — J Z z) 

j 2 = 

li 

4J±± — 4J± + 4a/2J 2 -j- + J zz ) 

h = 


'-4J±± - 2 J± + 2\[23 z ± - J zz ^j 

j 4 = 

h 

^23±± - 2 J± - 2V23 Z ± ~ Jzz) 
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FIG. 7. (Color online) (a) Experimental spectra obtained in reference [8] at H = 2 T and T = 0.05 mK. (b,c,d) Exam¬ 
ples of simulations, performed in the same temperature and field conditions, giving a good agreement with the data : (b) 
(J i , J 2 , J 35 J 4 ) = (-0.09,-0.19,-0.25,0.005); (c), (Ji, J 2 , J 3 , J 4 ) = (0.04,-0.29,-0.26,0.024), and (d): (Ji, J 2 , J 3 , J 4 ) = 
(-0.02,-0.34,-0.29,0.036) meV. 


B. Overview of the phase diagram 

The zero-temperature phase diagram for the nearest- 
neighbor pseudo spin-1/2 Hamiltonian given by Eq. © 
has been studied in Ref. [24]and[25l Owing to the XXZ 
nature of the Yb 3+ magnetic anisotropy, and depending 
on the sign and amplitude of the exchange parameters, 
four different magnetic phases are obtained in the phase 
diagram: a canted (or splayed) ferromagnetic phase (la¬ 
beled FM), the antiferromagnetic Palmer-Chalker phase 
which results from the so-called Y j representational, and 
the basis states ^2 and 7/3 of the Y 5 represent at iorP^, 
whose combination form the classically continuously de¬ 
generated manifold. Figure [10| shows a sketch of the 7 / 2 , 
'ips and FM spin configurations on one tetrahedron. 

As quoted previously, the parameters of the Hamilto¬ 
nian determined for Yb 2 Ti 207 in prior studies place it 
in the ferromagnetic phase. The aim of the next para¬ 


graphs is to test the unicity of this set of parameters, 
and further explore the consequences of the location of 
Yb 2 Ti 207 in the phase diagram. 


C. Exchange couplings in Yb 2 Ti 2 C >7 

As already reported in Ref. [ 8 |and 10, a polarized state 
sets in when a magnetic field H > H^ =0 05K = 0.5 T 
is applied along [1 I 0]. The spin dynamics consists in 
spin-wave excitations propagating in different directions 
of reciprocal space, as measured by inelastic neutron scat¬ 
tering for H = 2 and 5 T®. The data collected in these 
studies along directions (h h h), (1 1 £), (2 2 1), (h h 0 ) 
and (ft, ft, 1) are shown in Figure [T^a). Four modes can be 
distinguished in the experimental spectra, one of them 
being almost flat throughout the Brillouin Zone into the 
(ft ft 0 ) — (0 0 i) plane. 






































































To determine the exchange couplings, we sol ve the gen¬ 
eral Hamiltonian ([!]) in the RPA approximation 3 ^^, and 
perform exhaustive calculations of the spin-spin correla¬ 
tion function S(Q,uj). The quality of the fit relies on a 
good matching of the location of the maximum neutron 
intensity in (h h h), (1 1 t), (2 2 t\ (h h 0) and (h h 1) 
directions, as well as of the overall spectral weight 31 ! The 
calculated mean-field magnetic moment is checked to re¬ 
produce the experimental one 15 ! at T = 0.05 K for H = 2 
and 5 T. Finally JC a ,b,c,4 are transformed into Ji, 2 , 3 , 4 - 


Reduction of the parameter space by fitting the field-induced 
conventional spin-wave excitations 

In contrast with prior reports, we find that this pro¬ 
cedure does not lead to a unique solution. Representa¬ 
tive spectra are shown in Figure [7] (b,c,d) for different 
relevant sets. They all reproduce most of the experimen¬ 
tal features, yet slight differences, mainly in the spectral 
weight distribution remain. For instance, while the spec¬ 
tral weight of the two higher branches in the (lit) direc¬ 
tion is minimum at Q = (1,1,0), the inverse is observed 
in the simulations. More generally, the spectral weight 
distribution of the modes is never fully reproduced, al¬ 
though the agreement seems slightly better in the inter¬ 
mediate parameter range (e.g. Figure [7] (c)). 

The sets of parameters giving such good agreement are 
reported in Figure [8j Clearly, there is a linear relation be¬ 
tween the Ji s, so that the 4 parameters space is now re¬ 
duced to a one dimensional space. Colored areas delimit 
the uncertainties. Note also that, since both magnetic 
fields stabilize the same ferromagnetic ground-state, fit¬ 
ting the H = 2 and 5 T data provides similar constraints 
in the parameter space, resulting in the same linear re¬ 
lation between coupling constants. For this reason, the 
parameters obtained from the H — 2 and 5 T spectra 
are not distinguished and represented together in Fig¬ 
ure [8] It must be emphasized that the set of parameters 
reported in Ref. [8] is consistent with our determination. 

It should be stressed also that a number of parameters 
do not lead, in zero field, to the expected ferromagnetic 
phase (labeled “FM” in Figure [8]) but to the ^3 antiferro¬ 
magnetic ground-state. The phase boundary separating 
the 2 ps and the FM phases at T = 0, while exploring 
the one dimensional parameter space described above, is 
depicted in Figure [8] by a vertical dashed line and corre¬ 
sponds to Ji C = (—0.023, —0.326, —0.282,0.026) meV. 


Rods of diffuse scattering as a further constraint on the 
coupling parameters 

To further constrain the dimension of parameter space, 
other experimental results have to be taken into account. 
To this end, 5(Q,u; = 0) was calculated using a combi¬ 
nation of Monte Carlo and spin dynamics simulations 
at finite temperature, aiming at reproducing the rods of 



J 2 (meV) 

FIG. 8. (Color online) Exchange interaction sets determined 
based on the field-induced spin-waves at H = 2 and 5 T. 
Exchange interactions reported in Ref. [8]and[22jare indicated. 
Those reported in Ref. 9 and [16] are out of the parameter 
range of the figure, and are instead indicated in the Fig. [14] 
of Appendix [B] 


diffuse scattering. This numerical method is briefly de¬ 
scribed in Appendix [C] and more details may be found in 
Ref. |32] and [33] Calculations are based on the effective 
Hamiltonian © and were performed in the short-range 
correlated regime at T = 0.4 K, just above the critical 
temperature, for different coupling sets within the one 
dimensional parameter space determined above. Those 
numerical simulations (see Figure [9]) show that the rods 
of scattering around (2,2,2) and (2,2,0) are simultane¬ 
ously observed when approaching the FM/AFM phase 
boundary (in agreement with Ref. 25). This confirms 
the observation made in section [II B 1| of coexisting FM 
and AFM short range correlations. The cuts in direction 
(h h 2 — h) shown in Figure[3](b) and (c) further show that 
a good qualitative agreement between the experimental 
data (filled circles) and the spin dynamics simulations 
(empty circles) may only be achieved while lying very 
close to the FM/AFM phase boundary. 

Combining the different results, the 
best agreement is obtained for J^° = 


(-0.03(2), -0.32(1), -0.28(3), 0.02(2)) meV (grey 
area in Figure [8] (d)). The associated Curie-Weiss 
temperature is Ocw = 0.53(1) K, in excellent agreement 
with the e xperi mental value Ocw — 0.5 K determined 
in section |II A (see also Ref. therein). It is worth 
noting that this set of values locates Yb 2 Ti 207 closer 
to the AFM/FM phase boundary than the one obtained 
in Ref. [8] and [22] It is also quite different from that 
obtained by Thompson et aP as well as that reported 
by Chang et aPP. Actually, the latter studies yield spin 
wave spectra in strong disagreement with experiment, 
while our set reconciles both the diffuse scattering and 
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(hhO) (r.I.u.) 

FIG. 9. (Color online) Elastic scattering function S(Q,uj = 0) 
in the (h, h , l) plane obtained from spin dynamics simulations 
for J 2 = -0.26, -0.29, -0.31, -0.319, -0.326, -0.345 meV at 
T = 0.4 K. 


inelastic neutron data. 


(a) (b) (c) 



J 2 (meV) 

FIG. 10. (Color online) Sketches of the spin configurations on 
one tetrahedron in the (a) t/> 2 , (b) ^3 and (c) FM phases, (d) 
Phase diagram calculated as a function of temperature and 
of exchange coupling varying along the optimal one dimen¬ 
sional parameter space determined from the spin-wave fitting 
procedure. The black line roughly indicates the range of pa¬ 
rameters for which the transition toward the ferromagnetic 
order is first order. 




D. Zero field Ground-state. Order of the transition 

We now examine in more detail the consequences of 
the above results on the zero field ground state. To this 
end, the phase diagram within the ID parameter space 
discussed in the previous section is calculated based on 
Monte Carlo simulations, for different lattice size up to 
L = 16, N sp = 65536 (see Figure [f0|). Since the critical 
temperatures does not seem to evolve much with the lat¬ 
tice size for L > 12 (not shown), the phase boundaries 
have been roughly determined by locating the maximum 
of the specific heat for L = 12^. For the optimal val¬ 
ues Ji° close to the boundary located at J^ c , the specific 
heat and magnetic susceptibility show, upon cooling, a 
succession of several transitions. Two second order tran¬ 
sitions first occurs at about T = 0.275(3) and 0.174(3) K 
towards the ^2 and ^3 AFM states respectively. This is 
followed by an abrupt transition at T c = 0.15(1) K to¬ 
wards a FM state (sketches of the magnetic configuration 
are shown in Figure [To] (a,b,c)). The first order charac¬ 
ter of this second transition is demonstrated by a finite 
latent energy at the critical temperature. The latent en¬ 
ergy decreases while getting deeper into the FM state, 
such that the first order transition smoothly transforms 
into a second order one. This behavior while approach¬ 
ing the phase boundary is consistent with the study by 
Yan et alPP. 


Both the reentrant behavior of the FM phase at low 
temperatures and the evolution of the order of the tran¬ 
sition on approaching the phase boundary might result 
from thermal order by disorder: although the ferromag¬ 
netic state has a lower internal energy for J 2 > Jf, the 
^ 2/^3 phases may be favored in a small temperature 
range because of their strong antiferromagnetic fluctu¬ 
ations; the long-range ferromagnetic order is then stabi¬ 
lized at lower temperatures where thermal fluctuations 
are reduced. 


E. Spin dynamics simulations in zero field 

We finally turn to the calculations of the classical spin 
dynamics aiming at a final discussion of the experimental 
results. 


RPA approximation 

S(Q,uj) was first calculated in the RPA approxima¬ 
tion. Spectra have been averaged over equi-populated 
domains. They are presented in Figure [lljc) for differ¬ 
ent high-symmetry directions in reciprocal space. For 
comparison, experimental data in the same directions 
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FIG. 11. (Color online) S(Q,uj) along high symmetry direc¬ 
tions : (a) neutron data (this work) taken at T — 0.05 K; 

(b) calculated spectra obtained from classical spin dynamics 
simulations in the short-range correlated regime (T = 0.4 K); 

(c) RPA spectra obtained in the FM phase. 


are shown in panel (a). Calculations have also been 
performed for different exchange interactions around 
the FM/AFM phase boundary, along the lD-parameter 
space determined in the previous section. While the ex¬ 
citations are “gapless” in the AFM phase (we did not 
consider order by disorder or crystal field effects at this 
point, which could however be at the origin of a small spin 
gap), an anisotropy gap opens in the FM region because 
of a sizeable tilt of the spins from their easy-plane (see 
Figure [lljc)). Obviously, at this level of approximation, 
the model fails to reproduce the experimental data. Ac¬ 
tually, the calculated spectra consist in conventional spin 
wave excitations, which cannot capture the flat line shape 
of the experimental spectra. A more striking discrep¬ 
ancy especially affects the energy scale, which is about 
twice larger experimentally than in the simulations at 
most wave-vectors. In the simulations, the highest spin- 
wave branch is observed at Tiuj = 0.6 meV, a value which 
does not evolve much by varying the exchange parame¬ 
ters. Experimentally, the spectral weight reaches ener¬ 
gies Tiuj ~ 1 — 1.5 meV depending on the wave-vector 
(see Figure [4]). Note that such a difference cannot be ex¬ 
plained by the experimental resolution, which is around 
A E ~ 0.095 meV for kf = 1.15 A and an energy transfer 
hoo = 1 meV. 


Finite temperature spin dynamics simulations in the 
short-range correlated regime 

To go a step further, the classical spin dynamics was 
solved numerically by taking into account non-linear ef¬ 
fects associated with thermal fluctuations (see Appendix 
E2 and Ref. M and [33] for technical details about the 
method). Such spin dynamics simulations of S(Q,u) 
have been performed at different temperatures from T = 
0.1 to 0.9 K, on the basis of the effective Hamiltonian © 
and optimal exchange parameters Ji° proposed in the 
previous section. Intensity maps of S(Q,u) are shown 
Figure [TT] (b) for T = 0.4 K in high symmetry directions. 
The simulations succeed in reproducing the overall shape 
of the excitation spectrum observed experimentally: (i) 
in direction (h h 2 ), the excitation spectrum goes from 
quasi-elastic approaching Q = ( 2 , 2 , 2 ) to flat-toped ex¬ 
citations at ( 0 , 0 , 2 ); (ii) the excitations are mainly flat 
or weakly inelastic along direction (0 0 £); (Hi) in di¬ 
rection ( h h 2 — /i), a strong quasi-elastic signal over¬ 
whelm the spectra between Q = (1,1,1) and (2,2,0). 
More generally, the computed S( Q,o;) displays a quasi¬ 
elastic (resp. flat-toped/weakly inelastic) line-shape close 
to (resp. away from) the rods of scattering, as observed 
experimentally. This is illustrated Figure [l 2 j displaying 
the calculated scattering function for different temper¬ 
atures T = 0.3, 0.4, 0.5, 0.6 and 0.9 K, and wave-vectors 
Q = (1,1, 1 ), ( 2 , 2 , 0), ( 1 / 2 ,1/2, 3/2), and (0, 0, 2 ). While 
the two first positions are dominated by a strong quasi¬ 
elastic contribution, the latter are mostly inelastic. The 
experimental data measured at the same wave-vectors are 
shown in the insets. Moreover, the calculated spectra do 
not evolve much with temperature below T = 0.6 K. 
Experimentally, this behavior is even more impressive 
since the spectra do not depend on temperature between 
T = 0.05 and 2 K. 

However, despite this success, the spin dynamics sim¬ 
ulations fail to reproduce the large energy range, which 
does not exceed hoo ~ 0.6 meV in the calculations, as 
also observed in RPA simulations. 


IV. DISCUSSION 

In summary, our results show that for the optimal val¬ 
ues of coupling constants determined to reproduce the 
field induced spin wave spectra, classical calculations cap¬ 
tures some of the magnetic properties of Yb 2 Ti 207 . It re¬ 
produces different experimental features, such as the non¬ 
trivial structure of the diffuse scattering which points out 
coexisting FM and AFM correlations above T c , as well 
as the transition towards a long-range ordered canted 
ferromagnetic state at T c . At lower temperature, ther¬ 
mal AFM fluctuations have then been shown to affect 
the stability of the FM order: approaching the boundary 
from the FM side, the order of the transition smoothly 
goes from second to first order; sufficiently close to the 
boundary, a long ranged antiferromagnetic order may 
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FIG. 12. (Color online) Calculated S(Q,cj) obtained from 
classical spin dynamics simulations at different temperatures 
T = 0.3, 0.4, 0.5, 0.6 and 0.9 K, and for relevant wave-vectors 
Q = (1,1,1), (2,2,0), (1/2,1/2,3/2), and (0,0,2). Insets 
show the neutron data at the same wave-vectors. 


even be favored in a significant temperature range above 
T c . This succession of phase transitions observed close 
to the phase boundary could clarify the peculiar shape 
of the magnetization curve as a function of temperature 
on a single crystal recently reported (see Figure 2(a) of 
Ref. |T 5 j): with decreasing temperature, it shows first a 
reversible bump around T = 0.18 K, followed by the ir¬ 
reversible step at T c = 0.16 K associated with the long 
range ferromagnetic order. 

However, the classical approach fails to explain the 
strongly reduced ordered moment observed experimen¬ 
tally below T c as well as the unconventional nature of the 
spin dynamics, not only at very low temperatures, but 
up to To ^ 2 K. Both features are obviously connected, 
since the part of the moment which is not condensed in 
the static component remains, by essence, fluctuating. 

To explain this unconventional spin dynamics, the cou¬ 
pling to other degrees of freedom cannot be ruled out. 
This phenomena has been shown to b e of utmost im¬ 
portance in Tb 2 Ti 207 for instancj- 46 * 47 *, but there is, to 
our knowledge, no experimental evidence in favor of such 
a coupling in Yb 2 Ti 2 07 . Moreover, the unconventional 
dynamics appears around 2 K which coincides with the 
apparition of short range correlations, which is a clue for 
a purely magnetic origin. 

It is thus tempting to invoke quantum effects as the 
missing ingredient in the proposed model. Indeed, the 
proximity with the AFM magnetic states is likely a 
source of enhanced quantum fluctuations. To account 


for the unconventional spin dynamics, two scenarii come 
to mind: (i) the fractionalization of magnons into decon- 
fined new spin 1/2 particles and (ii) the scattering by 
longitudinally polarized two-magnons excitations. 

The former case is very well known in one dimensional 
system^. Here the so-called spinons, which can be 
seen as domain walls, propagate freely and separate two 
”Neel” configurations 49 50 . In two and three dimensions, 
the situation is more complex as the interactions create 
an attractive potential between spinons that bind them to 
form conventional magnons—. The quasi-degeneracy of 
the FM and AFM phases in the relevant part of the phase 
diagram phase for Yb 2 Ti 207 could cancel these interac¬ 
tions and release the spinons. In this view, the contin¬ 
uum would be characterized by the absence of edge sin¬ 
gularity at its lower threshold, with the exception of the 
”on rod” Q positions, as revealed by the quasi-elastic like 
shape. Actually, classical examples of continuum demon¬ 
strate that both situations are possible: the classical De 
Cloizaux-Person spectrum in the case of the spin 1/2 
Heisenbe rg ch ain, shows an edge singularity at the lower 
threshold 48 * 52 *, the spectra decreasing roughly as 1 /cj (the 
exact exponent depends on the anisotropy^ 3 ). However, 
for “XXZ” chains, close to the Ising limit of a strong 
coupling along the z-axis, have a spectrum determined 
by gapped solitons, with no singularity at the edgeiJ^l 
Theoretical studies have shown that the essence of the 
singularity relies on the fact that many low-lying modes 
add upP3. In this picture, recent theoretical work argue 
that the spontaneous magnetization along (001) below 
the critical temperature should confines_the spinons by 
creating a string tension between them 57 . This should 
result in the discretization of the two-spinon continuum 
into multiple spin-wave branches, a feature which is how¬ 
ever not observed in the data, that especially show no 
significant evolution below T c . 

Alternatively, since a strongly reduced ordered mo¬ 
mentum is the sign for significant non-linear longitudinal 
excitations, the scattering by a multi-magnons contin¬ 
uum should be considered. In this case, the ground state 
is long range ordered, with propagating magnons as the 
elementary (transverse) quasi-particles. In the longitu¬ 
dinal channel, however, fluctuations along the length of 
the ordered moment may developP^. Experimental real¬ 
izations are however quite scarce and usually give rise to 
a kind of tail above the spin wave dispersion which is still 
very well defined. For this reason, we believe that this 
second scenario is less relevant. 

It should be stressed finally that the measurements 
performed at T = 4.5 K (see Figure [6]) show that the 
spectra tend to become more conventional when increas¬ 
ing the temperature. Further inelastic neutron scattering 
measurements should be performed as a function of tem¬ 
perature, to look for an eventual transition or cross-over 
between quantum and classical regimes, with spin dy¬ 
namics respectively governed by quantum and thermal 
fluctuations. 
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V. CONCLUSION 

In this work, we show that the short range correla¬ 
tions that develop below 2 K in Yb 2 Ti 207 come along 
with non-conventional excitations: no well defined spin 
waves are observed in both the short range correlated 
and long range ordered regimes; rather, excitations spec¬ 
tra are characterized by a very broad and almost flat 
dynamical response which extends up to 1 — 1.5 meV 
and does not evolve with temperature below 2 K, some¬ 
times coexisting with a quasi-elastic response depending 
on the wave-vector. Based on Monte carlo Spin Dy¬ 
namics and RPA calculations, we determine a new set 
exchange couplings that allows to reproduce both the 
diffuse scattering in the short range correlated regime 
and the spin wave spectra observed in the field polarized 
phase. This careful determination of the exchange ten¬ 
sors places Yb 2 Ti 2 07 in the close vicinity to a FM/AFM 
phase boundary, and, in any case, far from the canonical 
quantum spin liquid state. This location points out the 
possible role of quantum fluctuations arising from these 
two competing phases. We show further that conven¬ 
tional RPA fails to reproduce the experimental excitation 
spectra at low temperature. Spin dynamics simulations 
performed in the short range correlated regime repro¬ 
duce however some features of the excitation spectra but 
lead to an energy bandwidth which is twice smaller than 
the experimental observations. We speculate that quan¬ 
tum fluctuations between FM and AFM phases govern 
the spin dynamics in Yb 2 Ti 2 07 and especially that the 
observed spectra rather correspond to a continuum of 
deconfined spinons as expected in quantum spin liquids, 
than simple paramagnons. 
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Appendix A: Neutron data analysis 

To further analyze the continuum of excitations ob¬ 
served in neutron experiments, the data can be modeled 
with the following function : 

/(cj) = C + G r (cj) + [1 + ti(cj)] ^ ^ Ri{uS) — Ri(-uj) 

i= 0,..,6 

(Al) 

C is a constant background, G r (uj) a standard Gaus¬ 
sian profile of width A r to model the incoherent elas¬ 
tic scattering at uj = 0, 1 + n{uj) is the detailed bal¬ 
ance factor and Ri(cj) a series of flat spectral bands 
with a fixed width of 2A r = 170 /ieV, and centered at 


-^= 0 , 1 , 2 , 3 , 4 , 5,6 — 0.19,0.36,0.53,0.70,0.87,1,04 and 1, 21 
meV (see top of Figure ^3] for an example): 

Ri(cu) = Ai if \uj — Ei\ < A r and 0 elsewhere (A2) 


The different spectra were analyzed through this model 
and the relative weights Ai/ V • Aj are reported in Fig- 
function of Q for the high symmetry direc¬ 


ure 


13 


as a 


tions. Those weights sequence in decreasing order as Q 
approaches the rods positions, close to (2, 2, 2), (1,1,1), 
(3/2, 3/2,3/2) and (2,2,0) (see the orange bars), while 
the weights i = 0,1, 2,3 become roughly equal away from 
the rod. The spectra then appear ’’flat toped”. Along 
(00/) a weak maximum can be detected for the i — 2 
band (E = 0.53 meV). 

These inelastic data, and especially the flat energy de¬ 
pendence demonstrates that on- and off-rods is a general 
feature of the spin dynamics throughout the Brillouin 
zone. Furthermore, the measurements show that the dy¬ 
namical response extends up to 1 — 1.5 meV. 


Appendix B: Connection between models 

Model in terms of the full magnetic moment 

The widely accepted starting point is given by the fol¬ 
lowing Hamiltonian 

H = H CE F + H exc + H z . (Bl) 

H z = gjJ2i H ■ Ji , is the Zeeman term, with H the ap¬ 
plied magnetic field and Ji the magnetic moment at site i. 
Hcef is the crystal electric field (CEF) Hamiltonian and 
H exc = JT. JiKijJj is a bilinear coupling Hamiltonian, 

where the interaction tensor /Q j couples next neighbors 
magnetic moments J at sites i and j. By symmetry argu¬ 
ments, the 9 coupling constants of the 3x3 tensor /Qj are 
reduced to only Here, we assume an exchange tensor 
which is diagonal in the (a, b, c) frame linked with a R-R 
boncPl 

J A = E J i + KbW + £*<%&) J j 

H,v=x,y,z 

1C 4 . y/2 bij . ( Ji X J j ) 

Alternatively, we can use 1C 1 , 2 , 3,4 defined by the simple 
relations: 

JC a + JC c 
/Cl = ^~ 

JC 2 = IC b 

„ )Ca-lC c 


Model in terms of a pseudo spin 1/2 

Since the energy gap between the CEF ground-state 
and the first excited levels are order of magnitudes larger 
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FIG. 13. (Color online) Top: example of raw data fitted ac¬ 
cording to the series of spectral bands Ri(cj). Bottom: Rel¬ 
ative weights Ai / JY Aj of the different bands for the data 
recorded at base temperature along (h h 2 ), (0 0 Q, (1 It) 
,(/i h 2 — h) and {h h 3 — h). Ai is defined via the spec¬ 
tral function (1 -\-n(cj)] JT =0 6 Ri(u) — Ri(-uj), where Ri(uj) 
is constant, equal to Ai in a given energy range centered at 
Ei= o,i, 2 , 3 , 4 , 5,6 = 0.19, 0.36,0.53,0.70, 0.87,1.04 and 1.21 meV, 
with a fixed width of 2A r = 170 geV. 


Full magnetic moment 
(K) 

Pseudo spin 
(meV) 

Pseudo spin 
(meV) 

JC a = -0.48 

J±± = 0.04 

Ji = -0.03 

K b = 0.23 

J± = 0.085 

J 2 = -0.32 

JC c = -0.67 

J z ± = -0.15 

J 3 = -0.28 

/C 4 = 0.02 

Szz = 0.07 

J 4 = 0.02 


TABLE I. Optimal sets of exchange in Yb 2 Ti 2 07 


The connection between JC and J is realized using the 
anisotropic Lande factor g±, z 30 . Defining Ax,^ = 


J« = 

J± = -A 


K a - 2/C c - 4/C 4 
3 

2 2/C a - 3/C b - K c + 4/C 4 

L i o 


Aj^ A z 

J±zh ^X 


tea ~\~ JC c — IC4 
3V2 

2 2/C a + 3JCb — JC c + 4 /C 4 
12 


Alternatively, a set of effective parameters Ji 2 3 4 can be 

usecPS 


Ji = - ^2J±± + 4J± + 2\/2Jz± — J zz^j 
J 2 = ^ ^4J±± — 4J± + 4\/2J z ± + J zz^j 
= - ^-4J ± ± - 2J ± + 2V%J Z ± - J zz'j 
Ja = - ^2J±± - 2Jx - 2\[23 z ± - 

Figure [14] provides the different sets of parameters and 
correspondence between these 4 conventions able to cap¬ 
ture the spin wave spectra measured in the field polarized 
phase ofYb 2 Ti 2 0 7 . The correspondence between the op¬ 
timal sets, obtained by taking into account the fit of the 
zero field diffuse scattering, is finally given in Table [TJ 


than the exchange interactions and the Zeeman term, it 
is possible to define effective spin 1/2 operators, denoted 
by S$, by projecting the full moment J i onto the CEF 
ground-state doublet. The effective spin 1/2 Hamiltonian 
reads: 

^eff = (B2) 

ij 

A popular convention consists in using (J±±, J±, J z ±,J Z z) 
defined as: 

H eS = Y J~S?S| - J ± (S+S- + S-S+) 

+J±± + 7**jSj Sj ) 

+J z± [Sf + QSj ) + * ++ j] 

where 7 ij,(ij are numbers defined in Ref ill 23, and 
M Note that ” sanserif’ notations refer to local bases. 


Appendix C: Classical spin dynamics simulations 

We consider the Heisenberg model H = Yh{ij) JijSj, 

where the summation is limited to nearest neighbors, Jij 
is the interaction tensor, and |S/ = 1/2 are classical 
pseudo-spins located at the pyrochlore sites. Our inter¬ 
est lies in the time evolution of the spin-pair correlations 
emerging in such a model. It is convenient to probe such 
dynamical correlations by calculating the dynamical scat¬ 
tering function S(Q, cu), which can be done by combining 
Monte Carlo and spin dynamics methods. 

The classical dynamics of the pseudo-spins is described 
by the non-linear Bloch equations 

f = s,x^4s, (ci) 

j 

where sites j are the nearest neighbors of i. These equa¬ 
tions of motions were numerically integrated using an 
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FIG. 14. (Color online) Correspondence between the different 
exchange coupling sets. The green, red and blue shaded areas 
correspond to the range of parameters giving a good agree¬ 
ment with the field induced spin wave spectra (see section 

IIIC). 


8th-order explicit Runge-Kutta (RK) method with an 
adaptive step-size control, offering an excellent compro¬ 
mise between accuracy and computation time. 


The initial spin configurations used for the numer¬ 
ical integration are generated at each temperature by 
a hybrid Monte Carlo method using a single-spin-flip 
Metropolis algorithm 32 33 . A thousand spin configura¬ 
tions are used at each temperature to evaluate the ensem¬ 
ble average in S(Q,uj) while the number of Monte Carlo 
steps needed for decorrelation is adapted in such a way 
that the stochastic correlation between spin configura¬ 
tions is lower than 0.1. The numerical results have been 
obtained for different lattice sizes ranging from L = 8 
(Ns = 8192) to L = 16 (Ns = 65536) with periodic 
boundary conditions. 


The classical Monte Carlo approach described above 
was also used to derive thermodynamics quantities such 
as the specific heat and the magnetic susceptibility, which 
allowed us to determine the phase diagram presented in 
Figure [lOj 
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